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In a recent work Zhang and coworkers (PRB 82, 075111 (2010)) studied the Zeeman sphtting of 
the Kondo resonance for the single impurity Anderson model in a finite magnetic field B with the 
numerical renormalization group (NRG) method. They report a discrepancy between the position 
of the Kondo resonance in the total spectral function and the position in the spin resolved spectral 
function at large magnetic fields. Additionally, the position of the Kondo maximum exceeded the 
Zeeman energy for B/Tk ^5 — 10, where Tk is the low energy Kondo scale of the model. In this 
comment we argue that both these findings are a result of the specific choice of NRG parameter 
values. However, we reproduce the crossover in the splitting from Kondo-like behavior to a non- 
universal splitting larger than the Zeeman energy, but this crossover occurs at much larger fields of 
the order of the charge scale. 



In a recent workQ Zhang and coworkers studied the 
Zeeman splitting of the Kondo resonance for the single 
impurity Anderson model (SIAM) in a finite magnetic 
field B. They employed the numerical renormalization 
group (NRG) to obtain spectral functions and discussed 
the position of the Kondo resonance in the spin-resolved 
and total spectral function, Pcr(w) and p{u!) = ^^^.^^(w), 
respectively. Two important points made in the work 
are: 

(i) With increasing magnetic field B the position 5-^ of 
the Kondo resonance in the total spectral function 
p{co) does not approach its position A-j- in the spin 
resolved spectral function, but instead > 5^. 

(ii) The positions S-^ and A-^ exceed the Zeeman energy, 
i.e. 5i~,A| > B for B /Tk > 5 - 10, where Tk is 
the low energy Kondo scale of the model {g = 2, 
fJ-B = ks = fi ^ 1). 

We show in the following, that both these conclusions 
cannot be drawn from the data presented in the publica- 
tion, since they originate from the specific choice of NRG 
parameters. 

We can indeed reproduce the published data[H for the 
SIAM at T = with [/ = -2e = 1 and T = 0.16, which 
are depicted as set ^1 in Fig. [1] For this set the dis- 
cretization parameter is A = 2.5, only a small number 
of states Ns = 150 are retained in each NRG iteration, 
a large broadening of a = 0.8 is used, = 20 dif- 
ferent discretizations are averaged (z-averaging), and an 
additional shift Q 7 = a/4 is employed in the broadening 
procedme. In the following, however, we demonstrate the 
strong sensitivity of peak maxima on the NRG param- 
eters. While set ^2 differs from set #1 only by setting 
7 = 0, for set ^^3, a discretization parameter A = 2 is 
used, Ns = 1800 states are kept, a = 0.075, and = 12 
different conduction band discretizations are averaged. 
All data are calculated within the complete Fock-space 
algorithm 3 also employed by Zhang et al. [l[ which co- 
incides with full density-matrix approachQ at T = 0. 

The spin-resolved spectral function p^{uj) in a finite 
magnetic field B = 0.003273 « 2.2Tk is shown in 
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FIG. 1: (Golor onUne) (a) Spectral function of the SIAM 
in a finite magnetic field for three different NRG parameters 
(see main text). The inset displays a close-up of the region 
around the Fermi energy and the arrows indicate the position 
of the maximum, (b) Maximum position of the Zeeman-split 
Kondo resonance in the spin-resolved (A^) and total spectral 
function (5^) as function of magnetic field, (c) as function 
of the magnetic field. The inset shows for set #3 the position 
normalized with the magnetic field. Ai^/B = 1 corresponds 
to the Zeeman splitting. The upper scale on panels (b) and 
(c) measures the field in units of the Kondo temperature Tk- 



Fig. [Ija) for these three different parameter sets. While 
the results share the same qualitative features, they still 
differ considerably. The differences do not only affect 
high energy features like the Hubbard satellites, but also 
the low energy Kondo resonance (see inset). The width 
of the latter is much narrower for set ^3 and, most im- 
portant for the present discussion, the positions of the 
maximum differ (indicated by the arrows in the inset). 
In panel (b) the positions of the Kondo maxima in total 
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and spin-resolved spectral functions are compared. For 
set #3, and very quickly approach each other and 
are indistinguishable for B > 0.004 w 3Tx (see upper 
scale in the plots). In contrast, the curves for set #1 
do not approach each other as observed in Ref. [l]. Panel 
(c) displays the maximum of the Kondo resonance A-|- in 
p^{uj) as function of the applied magnetic field B. The 
curve of set #3 increases almost linearly for all fields. It 
is slightly below the Zeeman splitting and exceeds it only 
for very large fields B > 0.2 = 1.25T « ISOT^ (see in- 
set). In contrast, the curves for sets #1 and #2 increase 
much faster and exceed the Zeeman-splitting already for 
fields B > 0.02 w 15Tk. 

The discretization of the conduction bandQ within 
NRG represents a severe approximation and recovering 
the correct continuum limit (e.g. A — > 1+ and Ns — > oo) 
is conceptually far from trivial. The broadening proce- 
dures adopted in the standard approaches 0, [H| offer no 
concise method, how to choose specific values for the pa- 
rameters involved. But as these parameters are artificial 
and remain A dependent, the physical phenomenon of 
interest must not depend on the actual value of these pa- 
rameters. In Fig. [T] we demonstrate that the results of 
Zhang et al. [ij are non-universal since they are strongly 
dependent on the broadening parameters. Even chang- 
ing only the artificial shift from 7 = a/4 (as in set ^1) 
to 7 = (as in set ^2) considerably changes the peak 
positions. In contrast, changing the same parameter for 
set #3 (within reasonable bounds, e.g. < 7 < a) does 
not change the spectral functions. This also holds for 
the other parameters, where the artificial broadening a 
is most critical. Increasing a reduces the critical field, 
where the A-|- exceeds the Zeeman splitting. Generally, 
the accurate extraction of Zceman-split peak positions 
from NRG spectral functions is a very delicate task as 



these always occur at finite energies, where the discretiza- 
tion errors are in principle not negligible. 

We therefore argue that the two results of Zhang et 
al. [l| stated above are inconclusive and have been ob- 
tained by a particular choice of NRG parameters. Since 
there docs not exist a "correct" choice for these parame- 
ters, physical conclusions must only be drawn from NRG 
spectra which arc robust to changes of NRG parameters 
or error bars should be given. Neither requirement ap- 
pears to be fulfilled by Ref. [l[ as demonstrated above. 

However, we reproduce the crossover in the Zeeman 
splitting from Kondo-like behavior, Af/B < 1, to a non- 
universal splitting with A-|- /B > 1 [see inset in Fig [Ijc)] . 
But this crossover occurs at fields of the order of the 
charge scale, B ^T, that is for fields two orders of mag- 
nitude larger than reported in Ref. P . 
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In a recent work[lj Zhang and coworkers studied the 
Zecman splitting of the Kondo resonance for the single 
impurity Anderson model (SIAM) in a finite magnetic 
field B. They employed the numerical renormalization 
group (NRG) to obtain spectral functions and discussed 
the position of the Kondo resonance in the spin-resolved 
and total spectral function, Pai^) and = '^a- P<^('^)^ 
respectively. Two important points made in the work 
are: 

(i) With increasing magnetic field B the position 6f of 
the Kondo resonance in the total spectral function 
p{uj) does not approach its position A-|- in the spin 
resolved spectral function, but instead > S-f. 

(ii) The positions 6-^ and A-^ exceed the Zceman energy, 
i.e. (5t,At > B for B/Tk > 5 - 10, where Tk is 
the low energy Kondo scale of the model {g — 2, 
fJ-B = ks = h = 1). 

We show in the following, that both these conclusions 
cannot be drawn from the data presented in the publica- 
tion, since they originate in an improper choice of NRG 
parameters. 

We can indeed reproduce the published data[l| for the 
SIAM at T = with [/ = -2e = 1 and T = 0.16, which 
are depicted as set #1 in Fig. [T] For this set the dis- 
cretization parameter is A = 2.5, only a small number 
of states Ns = 150 are retained in each NRG iteration, 
a large broadening of a = 0.8 is used, Nz = 20 dif- 
ferent discretizations are averaged (z-averaging), and an 
additional shift Q 7 = a/4 is employed in the broadening 
procedure. In the following, however, we demonstrate the 
strong sensitivity of peak maxima on the NRG param- 
eters. While set #2 differs from set #1 only by setting 
7 = 0, for set #3, a discretization parameter A = 2 is 
used, Ns = 1800 states are kept, a = 0.075, and Nz = 12 
different conduction band discretizations are averaged. 
All data are calculated within the complete Fock-space 
algorithm Q also employed by Zhang et al. [l[ which co- 
incides with full density-matrix approach^] at T = 0. 

The spin-resolved spectral function Pf(w) in a finite 
magnetic field B = 0.003272 k. 2.2Tk is shown in 
Fig. [I^a) for these three different parameter sets. While 
the results share the same qualitative features, they still 
differ considerably. The differences do not only affect 
high energy features like the Hubbard satellites, but also 
the low energy Kondo resonance (see inset). The width 
of the latter is much narrower for set 7^3 and, most im- 
portantly for the present discussion, the positions of the 
maximum differ (indicated by the arrows in the inset). 
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FIG. 1: (Color online) (a) Spectral function of the SIAM 
in a finite magnetic field for three different NRG parameters 
(see main text). The inset displays a close-up of the region 
around the Fermi energy and the arrows indicate the position 
of the maximum, (b) Maximum position of the Zeeman-split 
Kondo resonance in the spin-resolved (A^) and total spectral 
function {5^) as function magnetic field, (c) A^ as function 
of the magnetic field. The inset shows for set #2 the position 
normalized with the magnetic field. Af/B = 1 corresponds 
to the Zeeman splitting. The upper scale on panels (b) and 
(c) measures the field in units of the Kondo temperature Tk ■ 



In panel (b) the positions of the Kondo maxima in total 
and spin-resolved spectral functions are compared. For 
set #3, 6^ and A-|- very quickly approach each other and 
are indistinguishable for B > 0.004 w 3Tk (see upper 
scale in the plots). In contrast, the curves for set #1 
do not approach each other as observed in Ref. [H Panel 
(c) displays the maximum of the Kondo resonance A-|- in 
p^{uj) as function of the applied magnetic field B. The 
curve of set #3 increases almost linearly for all fields. It 
is slightly below the Zeeman splitting and exceeds it only 
for very large fields B > 0.2 = 1.25F « 150Tk (see in- 
set). In contrast, the curves for sets #1 and #2 increase 
much faster and exceed the Zeeman-splitting already for 
fields B > 0.02 « 15Tk. 

The discretization of the conduction band 3 within 
NRG represents a severe approximation and recovering 
the correct continuum limit (e.g. A — s- 1+ and Ns — ?► oo) 
is conceptually far from trivial. The broadening proce- 
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dures adopted in the standard approaches [J, |5| offer no 
concise method, how to choose specific values for the pa- 
rameters involved. But as these parameters are artificial 
and remain A dependent, the physical phenomenon of 
interest must not depend on the actual value of these pa- 
rameters. In Fig. [T]we demonstrate that the results of 
Zhang et al. [l[ are non-universal since they are strongly 
dependent on the broadening parameters. Even chang- 
ing only the artificial shift from 7 = a/4 (as in set #1) 
to 7 = (as in set #2) considerably changes the peak 
positions. In contrast, changing the same parameter for 
set #3 (within reasonable bounds, e.g. < 7 < a) docs 
not change the spectral functions. This also holds for 
the other parameters, where the artificial broadening a 
is most critical. Increasing a reduces the critical field, 
where the exceeds the Zeeman splitting. Generally, 
the accurate extraction of Zeeman-split peak positions 
from NRG spectral functions is a very delicate task as 
these always occur at finite energies, where the discretiza- 
tion errors are in principle not negligible. 

We therefore argue that the two results of Zhang et 
al. [ij stated above are inconclusive and have been ob- 
tained by an improper choice of NRG parameters. How- 
ever, we reproduce the crossover in the Zeeman splitting 
from Kondo-like behavior, A-^/B < 1. to a non- universal 



splitting with A-|- /B > 1 [see inset in Fig [Ijc)] . But this 
crossover occurs at fields of the order of the charge scale, 
B ^ T, that is for fields two orders of magnitude larger 
than reported in Ref. [l|. 
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